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We investigate sum-rules applying to the Raman intensity in a strongly correlated system close 
to the Mott transition. Quite generally, it can be shown that, provided the frequency integration 
is performed up to a cutoff smaller than the upper Hubbard band, a sum-rule applies to the non- 
resonant Raman signal of a doped Mott insulator, resulting in an integrated intensity which is 
proportional to the doping level. We provide a detailed derivation of this sum-rule for the t-J 
model, for which the frequency cutoff can be taken to infinity and an unrestricted sum-rule applies. 
A quantitative analysis of the sum-rule is also presented for the d-wave superconducting phase of the 
t-J model, using slave boson methods. The case of the Hubbard model is studied in the framework 
of dynamical mean-field theory, with special attention to the cutoff dependence of the restricted 
sum-rule, and also to the intermediate coupling regime. The sum-rule investigated here is shown to 
be consistent with recent experimental data on cuprate superconductors, reporting measurements 
of Raman scattering intensities on an absolute scale. 



I. INTRODUCTION 

Restricted sum-rules (relating the partial integral of 
the intensity up to a cutoff with a correlation function) 
are useful tools for analyzing spectroscopic information in 
strongly correlated materials. For example the sum-rule 
for the optical conductivityhas received a lot of attention 
in the context of cupratesjj @, H, 0, @ . 

In contrast the sum-rules for the Raman scattering in- 
tensity have been less studied. These were first consid- 
ered in 00]. Recently it has been stressed that, for a 
doped Mott insulator, the non-resonant Raman inten- 
sity should be proportional to doping[8] 9], provided the 
frequency integration is carried up to a cutoff which is 
such that contributions from the upper Hubbard band 
are not included. This proved instrumental in analyzing 
recent experiments on cuprate superconductors [9], that 
revived the debate about the relationship between the 
pseudo-gap and the superconductive gap in the under- 
doped phase of these compounds [lol [III . fl2l. Il3j. 

The goal of this paper is to analyze the Raman sum- 
rules and their dependence on the strength of the inter- 
actions, doping, temperature and choice of upper cutoff. 

In section Hill we show that for the t-J model the right- 
hand side that enters the sum-rule for Raman scattering 
is proportional to doping. In this case, the frequency 
cutoff can be taken to infinity (since the upper Hub- 
bard band is absent due to the constraint of no double- 
occupancy), hence making a rigorous theoretical analysis 
easier. 

We then evaluate explicitly the Raman response func- 
tion of a doped Mott insulator by solving the Hubbard 
model using Dynamical mean-field theory. We explore 
the region in which the Raman intensity scales with dop- 
ing and how this region varies with the cutoff used in the 



sum-rule, and contrast those results with that of a cor- 
related material slightly below the Mott transition. In 
section [V] we compare the results of our calculation, for 
different choices of the upper cutoff, with experimental 
data and then in section fVTl we conclude with predictions 
for the temperature dependence of the integrated low- 
energy Raman intensity, using a slave-bosons treatment 
of the t-J model. 



II. RAMAN SUM-RULES 

Raman scattering, is a photon-in photon-out process 
happening when an external electromagnetic field is ap- 
plied on an system, and its non-resonant cross-section 
can be calculated from the Fermi golden rule, reading: 



fl(q,fi) = 27T^ 



exp(-/3ei) 



(1) 



\g(ki)gQL f ) J2 44 (f\T rs (q)\i)\ 2 8(e f - Ei Si), 



where (3 is the inverse temperature, g(q) — (he 2 /Vujq) 1 ^ 2 
where V is the volume, ui q and e are energy and polariza- 
tion vectors of the photons (i, / indicate initial and final 
states of the process) and 



5 (q) = 



k 



d 2 e k 

dk r dka 



k+q/2,cr Ck -q/ 2 . CT ' 



(2) 



is the stress operator tensor. Sij is the energy of the 
initial or final state of the system, SI is the transferred 
energy, creates an electron of momentum k, is the 
one-electron dispersion of the model under consideration 
and r, s are cartesian components. 



2 



Then, by the fluctuation-dissipation theorem the scat- 
tering intensity in one channel 



/ rg (q,f2) = ^ 



eM J £i) \(f\r rs (ci)\i)\ 2 6(s f -e i -Q) 



U 

xy,(g,n) 

1 - e-f jn ' 



(3) 



is related to the imaginary part of a response function, 
i.e. 



Xrs (q,f)=ieW([r-(q,i),r'- s (0,0)]), 



(4) 



which is a stress-stress correlation function of the unper- 
turbed system. 

We will consider in the following only the q = con- 
tribution, since the photon momentum is always much 
smaller than the Fermi momentum. 



What is remarkable in Raman scattering 14] is that by 
tuning the polarization of the incident and of the detected 
outgoing photons one can exploit selection rules to sort 
out different processes in the probed material. This is of 
particular utility in order to separate the response due to 
electronic excitations in different areas of the Brillouin 
zone, thus allowing to probe k-dependent properties of 
the material. 

This is, in particular, a key issue issue for the physics 
of cuprates, in which nodal and antinodal regions of the 
Brillouin zone are known to behave in distinctly different 
manners. It is remarkable that, despite being a q ~ 
probe (as is optics), Raman scattering can still address 
momentum-selective issues by exploiting the dependence 
of the Raman vertex on the polarization of the incident 
and scattered light. Especially for a material with perfect 
in-plane equivalence of the a- and b-axis, the B 2g and 
B\g geometry mainly probe, respectively, the nodal and 
antinodal regions. 

In the non-resonant response, by choosing the direction 
of the polarizations e*, e*, one isolates contributions by 
the different elements r rs (q) of the stress tensor. 

As an example let's consider a square lattice with near- 
est neighbor and next-nearest neighbor hoppings t and t' , 
that give rise to the dispersion in k-space 

£k = — 2t(cos k x + cos k y ) + At' cos k x cos k y . (5) 

Then in the B 2g scattering geometry e l ,e^ are perpen- 
dicular between each other and parallel to the x-y plane 
axes, and the Raman scattering operator is r xy (q). From 
([2]) one then sees that the sum over the momentum space 
is weighted by the vertex factor 



d 2 e k 



= At' sin k x sin k 



(0) 



thus sorting out the contribution along the nodal direc- 
tions. 

Analogously one can show how in the B\ g geometry, 
with e l ,e^ perpendicular, but oriented along the lattice 
diagonals, one selects the antinodal directions in k-space. 



In this work we are interested in discussing the sum- 
rules fulfilled by the Raman B 2g non-resonant response. 
It is well known [7Llia|. that any susceptibility of the gen- 
eral Kubo form @ fulfills the following sum-rule: 



dfifixk, O j(fi) = <[0,[ir J o]]> J (7) 



where O is the operator to the fluctuations of which the 
susceptibility is associated. 

Hence, for the B 2g raman scattering this sum-rule 
reads: 



n Jo 



(8) 



Eq. © has been used in Ref. Q to normalize Raman 
spectra in mercury-based cuprates. In the following we 
will show that the right hand side of this sum-rule is 
proportional to doping in the t-t'-J model [22] ]. 



III. THE B 2g RAMAN SUM-RULE IN THE 
t-t'-J MODEL AND ITS PROPORTIONALITY 
TO DOPING 

In this section we report the analysis of the direct eval- 
uation of the r.h.s. of the B 2g sum-rule eq. ([5]) in the 
t-t'-J model on the square lattice. 

This model is represented by the hamiltonian 



H 



V C i,<7 C j,<? 



H.c.) +y^ y 4Jki(Sk ■ Si - -n k rii) 



(9) 

acting on the Hilbert space where states with doubly oc- 
cupied sites have been removed. 

Here cj is the operator that creates an electron with 
spin a in a state localized on site i, and t%j the hopping 
integral between single-particle states at sites i and j, 
Sk is the electron spin at site k and AJki > is the 
antifcrromagnetic superexchange and is zero for k = I 
and nonzero at (but not necessarily restricted to) nearest 
neighbours. 

The hopping matrix t%j is chosen to have nonzero el- 
ements only for nearest (Uj = —t) and next-nearest 
(tij = t') neighbors. The onsitc element is the chemi- 
cal potential term that tunes the filling (ta = —fJ>)- The 
bare electronic dispersion is given by eq. 

We have to evaluate the doping dependence of the av- 
erage of the double commutator, eq. ([5]), in this model. 

To implement the non-double-occupancy constraint we 
use the formalism of the Hubbard Operators. Definitions 
and some algebra tricks for this formalism can be found 
in appendix [X] We thus obtain an exact expression of 
the r.h.s. of the sum-rule. 

Both the stress operator ^ and the Hamiltonian H = 
H t + Hj can be readily expressed with these operators 
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X, 



a/3- 



where the stress component r xy has elements 
9 2 e k 

— = 41 

k 



7 = E< 



(10) 



H J - ^2 2Jkl ( X a/3 X f3a - X aa X /3l3) C 11 ) 
kl,a/3 



(12) 



= At' cos (k • r ij ) sin k x sin fc y . 

(13) 



that are nonzero only when i and j are next-nearest 
neighbour, as can be deduced by this expression (in order 
to lighten the notation we will drop the superscript xy in 
from now on). 



The evaluation of the commutators leads, for the hop- 
ping part [23j]: 



[[t XV , H t ],T XV ]] = ^ i( T abtbc - T bc t a b)Tbd - (T ac t c d - T cd t ac )T bc ] [X^X^X^X^ - X^X^X^X^] 
abcd.cx^(3 

+ ( T abhc - n c t ab )T cd [XZoD^D^Xqp - X$ Q D%pD h 0U7 X$ ^ , (14) 

abcd,a.(3a 



r 



where we have defined (see also appendix [X]) the bosonic 
operator D l af3 = (X^ + SapX^). All the other operators 
in the expression are fermionic. 

In the following analysis we will show that the average 
value of this expression and of the other parts of the 
double commutator is ~ 8 at the leading order. 

Indeed, after having brought together and collapsed 
all the operators that refer to the same site into one 
(X® XQ a — X£ a , see appendix [A]), expressions contain- 
ing Xq or fermionic operators vanish exactly as 8 — > 0, 
because of the constraint. 

This is an exact statement, but it can also be seen 
more explicitly in approximate schemes like slave-bosons, 
as we will show here. Indeed one can evaluate the doping 
dependence of the average of the different terms in the 
double commutator by considering that each fermionic 
operator on a different site carries a renormalization fac- 



tor proportional to b < 
a factor \b 2 \ ~ 8: 



y/S, while an operator Xq carries 



XaB— fiafi/3, X l aQ — f} a bi, XQ a —b i fi, 



Xl Q = b]bi 



(15) 

One can then evaluate, using these relations, the ex- 
plicit dependence in doping of each term of the sum cor- 
responding to the r.h.s. of ((SJ with the prescription of 
collapsing first the operators living on same sites into 
one. 

We will see that the coefficients of all terms containing 
no Xq and less than two fermionic operators on different 
sites (that we call "dangerous terms") vanish and thus 
the dependence is oc 8 at the leading order. 

Let's analyse then [[t xv , H t ], r xy ]]. 



For the first of the two four-fermion contributions in 
(I14p these " dangerous terms" occur when in the sum a = 
c and b — d. For example the first product in this case 
contributes to the sum (over a 7^ fe, a =/= (3) with terms 
like: 



vd yb ya 



yd yb 



-fLfacfLfbf} (16) 



since in this case one can anticommute the Hubbard 
operators in order to collapse the ones that refer to the 
same sites (X^X^ — I° a , see appendix IA"|) and then 
evaluate the averages using the slave bosons operatorial 
equivalences ([13)) . 

This term is "dangerous" since its average is not ob- 
viously proportional to doping. But it doesn't actually 
contribute to the sum since its coefficient in this case 
vanishes. Indeed for a — c and b = d its coefficient be- 
comes (r ab t ba - T ba t ab )T bb - (T aa t ab - T ab t aa )T ba = 0, since 

The second four-fermion product instead is dangerous 
when a = d and b — c and again its coefficient {r ab t bb — 

T bb t ab )Tba — {l~abt ba ~ T b at a b)T b b = 0. 

The two-fermion terms are dangerous when a = d and 
the coefficient becomes {r a btb c — T bc t a b)T ca . For this to 
be nonzero a and c have to be next-nearest neighbours. 
Then it is easy to see that on the cubic lattice (and within 
the ranges previously defined for the matrices and ry) 
it is impossible to chose the site b such that the coefficient 
is nonzero. 

Thus we have shown that all contributions of 
[[t xv , Ht], T xy ]] with a less than linear dependence in dop- 
ing vanish. 
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For the magnetic part of the hamiltonian one has: 



liT*y,Hj],T*y}} 



E 

abcd,a.^f3 

E 

abcd.crct^p 



^T a bT c d(Jbc - Jac - Jbd + J a d)(X^ o X o0 X^ Q X Oa - X% X 0a Xp X Q p) 
2T ab T bd (Jbc - Jac) [ X aO D a0 X O(r X pa + X< kp X p\P aa X K<j + X oO D aa X QP X< f3a 



I yc yd y\b ya ya r^b -yd yc yc ya T^tb yd yd T-tb ya yc yc yd T~ib ya 1 

+ A -af3 J ^a0 L/ /3a Ji -0a ~ ^aO^cra^Oa^ PP ~ A /3^ A aO L/ TO A 0tr — A o<3 U ao A <3a^ PP ~ A pp A aQ^ aa^Oa} ■ 

I 



(17) 



The first of the four-fermion terms is never dangerous, 
since even when a = d and b = c, the product X^X^ = 

Xq and thus its average is proportional to doping. 

The second of the four-fermion terms is dangerous 
when a — b and c = d. But in these configurations the 
coefficient becomes 2T aa T cc (J ac - J ac - J ac + J ac ) = 0. 

The analysis of the two-fermion terms is a little more 
involved. They are dangerous when a = d. The product 
in this case becomes (in D l af3 = {X l af3 + S a pXQ ) we can 
drop the Xqq operators [24|): 



ya yb yc > yc yb ya 
A ao A <j(i^Pa T A Q/ 3 A cra A /3cr 

I ya yb yC , yC yb ya 

+ A <x/3 A acr A /3« + A a/3 A /3<T A cra 

ya yb yc yc yb ya 

~~ ^ao^cra^-pp ~ ^PP^aa^pa 

ya yb yc yc yb ya 

^oa^ao^pp ~ ^PP^aa^aa- 



We used the fact that operators at sites a and d always 
commute with operators at site b since the coefficient 
vanishes when a = b or d = b. 

We then have to exchange the indices a and b in the 
III, IV, VII and VIII term, in order to match the spin 
indices with the remaining four terms (the operators in 
a and b commute). We obtain: 



ya yb yC 

^aa^ap^-pa 

ya yb yc 
^-aa^aP^Pa 

ya yb yc 
^aa^aa^PP 

Xa yb yc 
a<j^oa A PP ' 



XC yb ya 
aP^aa A pa 

yc yb ya 
^ap^aa^-pa 

yc yb ya 
A Pp A cja A Pa 

Xc yb ya 
Pp^aa^pa 



since the renaming of a and b changes the sign of the 
coefficient T a bT~bd(Jbc — Jac)- Thus the terms cancel two 
by two (I and III, II and IV, V and VII, VI and VIII). 

We can then conclude that all the terms that con- 
tribute to ([[t xv , H],t xv ]}) are at least a 5 and that the 
sum rule integral ([8]) is proportional to doping. Note 
that this is true irrespectively of the nature of the pos- 
sible long-range ordering (e.g. both in the normal and 
superconducting states) . 



IV. DMFT CALCULATION IN THE HUBBARD 
MODEL 

In this section we report the explicit calculation of the 
Raman response function for the doped Hubbard model, 



nil, 



(18) 



at strong and intermediate coupling within Dynamical 
Mean-Field Theory fl6| (DMFT). One of our main moti- 
vations for studying the Hubbard model is to investigate 
the restricted sum-rule, up to a frequency cutoff il c , and 
study the dependence on the cutoff frequency. 

Single-site DMFT freezes spatial fluctuations while 
fully retaining the local dynamics. It thus allows to calcu- 
late dynamical response functions in the local self-energy 
approximation. [25j | 

Since the Raman vertex is odd under e.g. k x — > 
— k x , vertex corrections vanish in the local self-energy 
approximation [l6j and the stress-stress response function 
eq. (fj| at q = reduces to the simple bubble-diagram of 
dressed propagators, that reads (in the imaginary time 
formalism) : 



X(ifi) 



d 2 e k . 
8 ^ x dk x dk y ' 



2 G(k,iv)G{k,iis + iQ), (19) 



and analytically continued, 



X "(fl) = I deV{e) I dujA(e, u)A(e, uj+n) [f(uj) -f(uj + Q)] , 

(20) 



where 



v(«) = £( 



<9 2 e k 



(4t') 2 sin2 k * sin2 M( £ - £ k) (21) 



is the appropriate form factor obtained by summing over 
k the product of two B 2g geometry vertices eq. ©• It 
can be calculated once and for all and for the dispersion 
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(|5]) is a smooth function. A(e, u>) is the spectral function 
A(k,u) = -~ImG(k,u) 

7T 

1 ImS(w) 
_ 7T (w + Jl - e k - i?eS(w)) 2 + (im£(w)) 2 ' 

(22) 

for ek = e and is easily accessible in single site DMFT(26l|. 
through the knowledge of the local self-energy (27j. 

Hence in this approximation the non-resonant Raman 
response is readily calculated by convoluting two single 
particle spectral functions with the appropriate kernel 

We have solved the Hubbard model with DMFT at zero 
temperature and finite doping for U /W = 2.5 (W is the 
bandwidth), that can be seen as a rather strong coupling, 
and for an intermediate coupling, U/W = 1.35. 
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Figure 1: (Color online) Evolution of the B29 Raman response 
function for U /W = 2.5 with doping ranging from S = 0.01 
to 5 — 0.25. The low energy feature increase in weight for 
increasing doping, the high energy one decreases. 

Let's first focus on the strong coupling case. Indeed 
this correspond to a doped Mott insulator, the critical U 
at half-filling for the chosen density of states (see below) 
being U c /W ~ 1.5, and at U/W = 3 the Hubbard bands 
being split apart by a gap A ~ 1.5W [281 ]. 

In figure [1] we show how the response function cal- 
culated in DMFT shows two principal features. These 
are expected when A > W from the general features of 
A(e, uj) and the smooth form of V(e). 

At zero temperature the Raman response formula sim- 
plifies further in: 

X//(0) = [ deV(e) [ duiA(e, w) A(e, u + Q) (23) 
J J— a 

The integral will be nonzero when, for some e, the two 
spectral functions shifted by Q from one another have 



an overlap in the region [ — 0, 0] , i.e. between the two 
quasiparticle resonances. 

Hence at low f2, yll will show contributions from the 
overlap of the two functions within the same Hubbard 
band, the one in which the quasiparticle peak lies. 

At higher il > W, since the width of the Hubbard 
band is ~ W, if the gap is bigger than the bandwidth 
(A > W), the two spectral functions will have negligible 
overlap. Indeed the Hubbard band of one will fall onto 
the gap of the other, and the response will be nearly zero. 
Then for even bigger f2 the two spectral functions will be 
shifted enough so to have the upper Hubbard band of one 
overlapping the lower Hubbard band of the other, and a 
finite response is found, thus generating a second, well 
separated feature at these higher frequencies. 

If instead A < W, this separation of the Response will 
in general not occur. 

The evolution of this response function with doping is 
such that the lower energy feature grows with increasing 
doping, while the higher energy one decreases. This is 
expected, since upon doping the system dilutes and be- 
comes less and less correlated thus causing a transfer of 
spectral weight back from high to low frequencies. 

Physically the separation of the response function into 
two features means that the Coulomb repulsion is large 
enough to separate in energy the particle- hole excitations 
happening in the Hubbard band of empty sites from the 
ones involving doubly occupied sites. 

In the strong coupling limit and at low doping this 
model can be mapped on a t-t'-J model and this allows 
us to identify the low- frequency response with that of the 
t-J model, so that we can use it to confirm our analysis 
on the sum-rule integral. 

We thus consider the response under a cutoff fl c ~ W, 
that includes only the processes of the low-energy feature 
of the Raman Response function, namely 

w(n c )= <msi x "Hub(fy^ dnn x " u (n),tfn c >w. 

Jo Jo 

(24) 

In fig. [2]we plot the calculated W(£l c ) as a function of 
different cutoff frequencies fl c . The highest value, £l c = 
1.2W fully includes the low-energy feature of the Raman 
Response in our case, thus reproducing the t-t'-J-model 
sumrule. 

It appears that the linearity in doping of VF(f2 c ~ W) 
holds well up to dopings of about 10%, then higher order 
terms start contributing considerably. 

By reducing the cutoff frequency one can see that the 
integral is still linear in doping. This is somehow ex- 
pected, since the whole lower feature scales with dop- 
ing, and vanishes as the insulating phase is reached for 
5 = 0. The linearity region actually increases and extends 
up to 15% doping (i.e. the whole underdoped region in 
cuprates). 

The smaller extent of the linear region for the higher 
cutoffs is probably due to the proximity of the high- 
energy feature in the response function, that would in- 
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Figure 2: (Color online) Integrated spectral weight up to 
a cutoff Q. c = 1.2W, that fully includes the lower energy 
feature) of the first moment of the Raman response func- 
tion, plotted versus doping. As predicted for the t-J model 
W(8) oc 8. Integral up to lower cutoffs are also plotted. 



deed be completely absent if we were to take the actual 
U/W -> oo limit. 

But in fact considering a cutoff lower than fl c ~ W is 
actually directly relevant to the analysis of experimental 
data on cuprates, as we will see in section [V] 

Let's now consider the intermediate coupling case. 

We studied the model at U = 1.35W = 0.9U c2 and 
obtained the results plotted in the upper panel of fig. [3] 
One imediately sees that the two features of the strong 
coupling case have now merged as the separation of the 
two Hubbard bands is less than W. 

Nevertheless, the low energy part of x"(fy still scales 
with doping. 

In the lower panel of fig. [3] we plot the value of the 
integral W(Sl c ) for small cutoffs as a function of doping. 
Although the intercept is indeed nonzero, owing to the 
fact that a quasiparticle resonance is still present at half- 
filling in this case (that testifies the metallic behavior at 
half-filling that one obtains when U is smaller than the 
Mott transition critical value U C 2), this intercept is prac- 
tically negligible for the lower cutoffs, given the very little 
weight of the quasiparticle resonance very close to TJ C %. 
Thus the behaviour of the sum-rule integral for small cut- 
offs is still ~ 5 in the underdoped region for U ~ 0.9f/ C 2. 

Indeed by raising the cutoff to values closer to the 
bandwidth one eventually loses the linear behaviour, ow- 
ing to the presence of the low-frequency tail of the spec- 
trum of transitions between the two Hubbard bands. 
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Figure 3: (Color online) Upper panel: non-resonant B^g 
Raman response function for U/W — 1.35 < U C 2- Different 
curves are, from bottom to top, for increasing dopings 8 = 
0.01 -5- 0.25. The two separate features present at large U/D 
have now partially merged, but the curves still scale with 
doping at low frequency. Lower Panel: sum-rule integral for 
cutoffs fi c — 0.165W (lower curve) and Q c ~ 0.375W (upper 
curve). Inset: sum-rule integral for all the chosen cutoffs in 
Fig. |1 from bottom to top Q c /W = 0.165,0.375,0.65,1.2. 
An intercept is clearly non negligible for Q c — 0.651T, the 
behaviour is no more linear for O c = 1.2W of the order of the 
bandwidth. 



V. CONNECTION WITH EXPERIMENTS: 
CUTOFF FREQUENCY AND COUPLING 
STRENGTH 



Some experimental setups are believed to give 
measures that are accurate enough to extract reli- 
ably absolute scattering intensities. We have com- 
piled and anayzed experimental data from two such 



references [17| L8|] , which are shown in figure |4j The linear 



dependence of the absolute intensities in the underdoped 
region (at least up to ~ 10% doping) is clearly visible for 
small cutoffs, up roughly to 2000cm" 1 ~ 0.25eV. 



7 



Hence if the experimental cutoff is as low as il c < 
2000cm" 1 ~ 0.25eV Raman spectra can indeed be nor- 
malized using a linear scaling in doping, as was done in 
Rcf. [|. 
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Figure 4: (Color online) Sumrule integral W(il c ) calculated 
from the experimental data on (Yi- y Ca y )Ba2Cu30e+x in the 
normal phase, taken from refs [13] and where absolute in- 
tensities of Raman scattering have been measured. The plots 
are rescaled by normalizing to the square of the cutoff fre- 
quency. The sum-rule linearity holds clearly up to S ~ 0.1 for 
both lower cutoffs. The line is a least square fit of the £l c — 
1000cm" 1 datapoints, discarding the S = 0.01 point that cor- 
responds to the antiferromagnetic phase. The curves for these 
lower cutoffs extrapolate linearly to zero or to a very small 
intercept. The curve obtained for Q c = 4000cm" 1 - 0.5eV 
instead doesn't seem to extrapolate to zero. As discussed in 
the text, this may indicate either that cuprates are in an in- 
termediate coupling regime [l9|. or that additional magnetic 
contributions (e.g. two magnons) contribute to the Raman 
signal in that frequency range. 

By raising the cutoff around 4000cm" 1 ~ O.beV one 
finds a non vanishing intercept as is visible in fig (31 

It is also very interesting to notice the nonvanishing 
intercept found for higher cutoffs. This is consistent 



with the very recent analogous findings of Comanac et 
al. in Ref. [19J, where by analyzing optical conductivi- 
ties of a wide range of cuprates the mid-infrared spectral 
weight (f2 c = 0.8eV) is found to be non vanishing de- 
spite a vanishing Drude spectral weight (fi c = 0.2eV). 
To explain these experiments these authors propose to 
model the optical conductivity of cuprates with a Hub- 
bard model whose interaction strength is less than the 
single-site DMFT U c2 ~ 1.5W. This is the value at which 
a pure Mott transition takes place, without additional ef- 
fects due to magnetic correlations. 

Indeed the calculation we performed in section IIV1 at 
U = 1.35W < U C 2, reproduces the behavior found in 
the experiments we reviewed. As is clear in the lower 
panel of fig. [3] while for intermediate cutoffs an intercept 
is sizeable, for the lower cutoffs it is instead very small. 
Moreover the antiferromagnetic correlations neglected in 
this single-site DMFT approach are likely to depress the 
low- frequency data, thus further reducing the intercept. 

There are two possible interpretations of this obser- 
vation of a finite intercept for the higher values of the 
cutoff (keeping in mind however that the Raman data 
used here are taken at slightly different temperatures, 
and that more data are needed). 

The first one is that cuprates are actually in the regime 
of intermediate correlation strength, as recently advo- 
cated by Comanac et al.[19( on the basis of a similar 
observation from optics. Another possible interpretation 
of the nonvanishing intercept found for fl c — 4000cm -1 
are additional contributions to the Raman scattering of 
spin origin, such as two magnon processes. To settle this 
issue clearly further work is needed. 

VI. THE B 2g SUM-RULE IN THE 
SUPERCONDUCTING PHASE OF THE t — t' — J 
MODEL: SLAVE BOSONS 

In this paragraph we calculate explicitly the right-hand 
side of eq. ©, by averaging the calculated expressions 
for the commutators, eq. (fl4|) and (fT7| . in the Slave- 
bosons mean-field formalism of Ref. [2(| used in the early 
studies of High-Tc's, that is the simplest approximation 
that allows to access the d-wave superconducting phase 
of the t-t'-J model, beside the normal phase. 

We recall here briefly this approach. The hamiltonian 
of the t-t'-J model, eq. lives in the restricted Hilbert 
space in which all states with a doubly occupied site are 
projected out. This can be expressed by means of a mixed 
fermion-boson hamiltonian: 

h = E Mlwlfj<> + ' - v E fl^ 

i,j,<j i<y 

+ 4J«[(5, • Sj) \(1 b\bi)(l b}b 3 )} 

+ E A *(E/iU' + & &- 1 )' ( 25 ) 
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in which the boson operators b\ (hi) represent the cre- 
ation (destruction) of a hole on site i, and thus carry the 
charge degrees of freedom, while the pseudofermions j ia 
(fia) carry the spin information. This hamiltonian lives 
in the enlarged Hilbert space represented by all the pos- 
sible states of the bosonic and pseudofermionic fields. To 
complete the mapping a constraint has to be introduced, 
in order to restrict again this Hilbert space to the physi- 
cal one, by excluding all the non-physical states. In doing 
this, one can also implement the no-double-occupancy 
condition, using namely: 

E/il/- + ^ = 1 (26) 

This is expressed in eq. (I25[) by means of the Lagrange 
multipliers Aj, that enforce the condition on every site. 

One then performs an Hartree-Fock-Bogoliubov decou- 
pling and a static mean field (such that A; = A and 
the constraint is satisfied on the average, i.e. (bjbj) = 

1 — Y^aifiafia) = ^) on this Hamiltonian by introducing 
the following mean-field parameters: 

K = 3J(J2flfi+ x ,a)/2 

cr 

A x = —A y = A = 3J\ftffi+xl — filfi+xl)/2 

(btb j ) = \{b i )\* = l- y £(flf iir )=8 (27) 

a 

We choose the phase in which K is isotropic (i.e. K x = 
K y = K) while the superconducting order parameter has 
a d-wave symmetry (A z — —A y ). This has been shown 
to be the actual ground state at finite doping between all 
the possible symmetries in this kind of mean-field 20]. 

The obtained mean-field Hamiltonian reads: 

Hmf = £ Ck/L/to + (E A k/k T /-ki + H.c.), (28) 

kcr ka" 

where 

£k = — 2(St — K)(cos k x + cos k y ) + 48t' cos k x cosk y (29) 

Ak = 2A(cos k x — cos k y ) (30) 

Eqq. (EU, ©, (Hi and ^ are the mean-field equa- 
tions, that have to be iterated until a self-consistent so- 
lution is found. 

The hamiltonian (|28|) is of the standard Bogoliubov 
quadratic form, and thus it is trivially solvable (by means 
of a Boboliubov substitution) One can then calculate any 
fermionic average by means of a Wick factorization in 
terms of the quadratic averages 

(/i/^> = E cos ( k -^)(/k CT /^>, 

k 

ifiafjs) = E C0S ( k ' r ij){haf-ka)- (31) 
k 



This mean- field hamiltonian of the t-t'-J model has the 
phase diagram shown in Fig[5]for t = 10 J and for t' = 
(Kotliar-Liu result) and t' = 0.3t (the value generally 
used for High-Tc superconductors like BiSCCO). 

For each doping there are two phases, a low temper- 
ature (superconducting) phase with A ^ and a high- 
temperature (normal) phase with A = 0. The transition 
temperature is denoted Trvb and is a decreasing func- 
tion of doping. The bosons are treated here as fully con- 
densed, because {b\bi) = \ (bi)\ 2 = 8 at all temperatures. 
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Figure 5: (Color online) Doping-temperature phase diagram 
of the t — t' — J model within the Kotliar-Liu slave-bosons 
mean-field approximation for J = O.lt and t' = and 
t' — 0.3i. The actual superconducting region is below a con- 
densation temperature of the bosons, where phase coherence 
sets in. This temperature is an increasing function of doping, 
so that the superconducting region is dome-shaped [2C| . 

If we were to allow bosons decondensation (which re- 
quires to include the link variational parameter (bfbj) as 
well), another crossover temperature (zero at half-filling 
and rapidly increasing with doping) would delimit a low 
temperature phase with phase coherence from the dis- 
ordered phase of incoherent pairs. The actual supercon- 
ducting region has a dome- like structure centered around 
a region of " optimal" doping where the two characteristic 
temperatures cross. 

In this formalism it is then possible to calculate the 
r.h.s. of eq. ([5]) by taking each Hubbard operator prod- 
uct in the two parts of the commutators, eq. (|14|) and 
(fl"7)) and obtain the expression in terms of the f's and the 
6's using eqq. (fT5)) . Then one performs the Wick factor- 
ization of each product (i.e. following the example of eq. 
(HU) one obtains (fl a f aa fLM = (fLfa a )(f br7 fb*) - 

(fLfba) (fLfaa) - (flafL) (fLfbo) ) and by means of the 
quadratic averages eq. (|3 1 [) calculated with Hmf within 
the converged solution one obtains the final result. 

This result is plotted in FigO as a function of doping 
at different temperatures. 
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The results shows that in the superconducting phase 
at T < Tjivb the sum-rule integral is linear in doping 
until 30-40% doping, thus allowing the linear scaling of 
Raman spectra in doping for the full extent of the super- 
conducting phase of cuprates. 

Still one has to be aware of the many limitations of 
this simple mean-field. Beside the cited boson deconden- 
sation, that adds a second energy scale to the problem, 
less relevant to our analysis, the method lacks a second 
fermionic energy scale that is the new and very debated 
point put forward in Our analysis in terms of the 
Kotliar-Liu slave-bosons has to be taken thus as a qual- 
itative description of the nodal physics in the supercon- 
ducting phase. 
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Figure 6: (Color online) Right hand side of the sumrule evalu- 
ated in the t—t' — J model within the Kotliar-Liu slave-bosons 
mean-field approximation for different temperatures and pa- 
rameters J — 0.1t,t' — 0.3t. The bottom curve is calculated 
for f3 = J, entirely in the normal phase, and shows how the 
linearity holds pretty well also for the normal phase in the 
whole range of dopings. 

It is also interesting to plot the temperature depen- 
dence of the sumrule integral. 

In order to compare the temperature dependence in 
the superconductive phase at different dopings, in fig. [7| 
we plot ([[t xv , H},r xy ]) /S, owing to the doping linearity. 
The rescaled curves indeed lie all in the same range of 
values but show a different behavior in temperature de- 
pending on the doping value. While at low doping raising 
the temperature causes a decrease of the sumrule value, 
after 5 ~ 12% (that could be interpreted as "optimal 
doping" in this Slave Bosons mean-field ) the value in- 
stead increases with increasing temperatures. Then in 
general when the system enters the normal phase upon 
heating the integral value shows a cusp and then drops 
quickly. 

This behaviour has not yet been measured and it is a 
prediction of this slave boson mean-field treatment that 
would be very interesting to clarify in further experi- 
ments. 
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Figure 7: (Color online) Temperature dependence of the sum- 
rule in the Slave bosons formalism in the superconducting 
phase. For increasing temperature the integral decreases at 
low doping while it increases at higher doping. 



VII. CONCLUSIONS 

We have studied a restricted sum-rule for Raman 
scattering, which involves the first moment of the non- 
resonant Big response function up to a cutoff frequency: 

For the doped t-t'-J model, where we can take the 
cutoff O c = oo, we have shown that this integral is pro- 
portional to the doping level. This supports the experi- 
mental data analysis of Ref. @ . 

We have then calculated the Raman response function 
and studied its doping dependence for the normal phase 
of the Hubbard model within DMFT. 

We have first studied the system in the strong cou- 
pling case U = 2.5W, that corresponds to a doped Mott 
insulator. In this case the low- frequency response of this 
system can be identified with the one of the t-t'-J model, 
and we have confirmed its linearity in doping at low dop- 
ing. 

We have then studied the intermediate-coupling case 
U = 1.35W = 0.9U C 2, showing that in this case the sum- 
rule integral still shows a behavior ~ S for very low cut- 
offs, while the intercept is non-negligible for intermediate 
cutoffs. 

We have also shown a compilation of experimental data 
that are believed to correctly measure the absolute in- 
tensity of Raman scattering on YBCO. As in our DMFT 
study at U = 1.35W, the behaviour is ~ 6 at low cut- 
offs (f2 c < 0.125el/) but shows a non- negligible intercept 
already for fl c ~ 0.25eV. This can be interpreted in sup- 
port of the hypotesis that cuprates are "intermediately 
correlated" as recently put forward in Ref. 

Finally we have studied the superconductive phase of 
the t-t'-J model using the slave-boson method of Ref. 
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[20] |. We can calculate explicitly the right-hand side of the 
sum-rule in this approximation and we show its linearity 
in doping up to dopings of order of 30%, and a prediction 
on the temperature dependence. 

The combination of those results give a first indication 
on how the strong correlations affect integrated Raman 
intensities. Missing at this point, a systematic study of 
the effect of dynamical short range magnetic correlations, 
since the slave boson method captures these in a prim- 
itive fashion using static link expectation values. These 
limitations can be removed using cluster versions of dy- 
namical mean field theory[2l| and this is left for future 
studies. 
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Appendix A: SOME ALGEBRA OF THE 
HUBBARD OPERATORS USED IN THIS PAPER 

Hubbard operators on site i are defined as: 

X* af} = \aUl3\u |a),|/3) = |0>,|T>,U) (Al) 

and hence they are used to enforce the constraint of no 
onsite double occupancy, since they all project the | tl) 



state to zero. 

From their definitions it is easy to show that two Hub- 
bard operators acting on the same site observe the fol- 
lowing commutation (anticommutation) relations: 

X js}± = &Pl X a& ± 8a5X*p (A2) 

Instead two operators acting on different sites respect 
the canonical commutation or anticommutation relations 
depending on them being "fermionic" (X 0t7 or X a0 that 
add or remove a spin-1/2) or "bosonic" (X o, X aa > that 
either don't change the onsite spin or they add or remove 
a spin 1), i.e. 

[x* B) X B ] = 0, [X* p , x F ) r 0, [X* B , X F ] = 0, V i j 

(A3) 

where " B" is for " bosonic" and " F" is for " fermionic" . 
Indeed these operators are not actual bosons or fermions 
since on-site they respect the (|A2|) . 

It is useful in our calculations to introduce the 
"bosonic" operator D % a ^ corresponding to the particu- 
lar case of (|A2[) for the onsite anticommutation of two 
fermionic operators: 

[Aq ! X l)f}] + = X a(3 + $al3 X 00 = D af3 (A4) 

and also its on-site commutation relations with the 
others Hubbard operators: 
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